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By deriving a generalized Harper's equation, we investigate the effects of temperature T and hole 
doping 8 on antiferromagnetically correlated electrons under the influence of an applied magnetic 
field. We obtain a phase diagram in a T-5 plane for staggered magnetization under the external 
field. We find that away from half filling the maximum values of staggered magnetization occur 
with rapid convergence to a finite temperature constant near half filling. We discover the reentrant 
behavior of the staggered magnetization in the presence of the magnetic field. 



PACS numbers: 75.50.Ee, 71.10.-w, 71.27.+a, 71.10.Fd 

and using the Landau gauge A = B(0, x, 0), we obtain, 
in the momentum space, 



Hubbard model is one of the most widely used model 
for studying the systems of correlated electrons. The 
original Harper's equation [^]-^). does not contain the 
effects of electron correlations and thus it takes into 
account only non-interacting electrons under external 
fields. To the best of our knowledge, there is no study of 
Harper's equation which takes into account such correla- 
tion effects. By considering an infinite square lattice, we 
derive a generalized Harper's equation which allows to 
investigate antiferromagnetic electron correlations under 
the applied magnetic field. Earlier, our study was lim- 
ited only to zero temperature and half filling [Q. Here 
we study how the applied magnetic field affects the sys- 
tem of antiferromagnetically correlated electrons at finite 
temperatures and away from half filling. 

The Hubbard model describing the two-dimensional 
system (square lattice) of electron correlations under an 
external magnetic field is given by 
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where t is the hopping integral; A, the electromagnetic 
vector potential; <fio — — , the elementary flux quan- 
tum; U, the on-site Coulomb repulsion energy, and /i, 
the chemical potential, (ij) stands for summation over 
nearest neighbor sites i and j. c| CT (c icr ) is the creation 
(annihilation) operator of an electron of spin a at site i, 
and ritf (n^i ), the number operator of an electron of spin 
up (down) at site i. 

We are interested in a uniform staggered magnetization 
(antiferromagnetic (AF) order), 

m = e iQri (ctfCj 
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where Q = (7r,7r), and = (i Xl i y ) with i x and i y be- 
ing integers including 0. The lattice spacing is set to 
be a unity. Introducing the mean field (Hartree-Fock) 
approximation, 
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where g= (2vr^,o) = (2*2,0 
flux quanta per plaquette. For simplicity we considered 
a uniform doping 8 = 1— {nA in the expression above. 
The first bracketed term in (H) represents hopping pro- 
cesses; the first term in the bracket is contributed from 
the nearest neighbor hopping in the x-direction and the 
last two terms in the bracket, from the nearest neighbor 
hopping in the y-direction. The external vector poten- 
tial A shifts the wave vector of electron by g during its 
hopping in the y-direction. The second term results from 
the antiferromagnetic electron correlations, and the third 
term shifts the total energy of the system as a result of 
hole doping. We do not consider the magnon excitation. 
Thus the exchange terms in (^|) above vanish due to the 
global SU(2) symmetry. 

The diagonalization of (^) yields the following general- 
ized Harper's equation which now incorporates electron 
correlations, 
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where by setting 


M n = -2 cos{k x + 


ng) and t 
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Here the matrix T is associated with electron hopping, 
involving phase modulation due to the external field, and 
U, with the antiferromagnetic electron correlations. 
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FIG. 1. Reduced Brillouin zone (shaded region) for the 
computation of energy dispersion relation with (a) half a flux 
quantum per plaquette, £ = |, and (b) no magnetic field, 
f = 0. Here the elementary flux quantum c/>o is set to be a 
unity. 

Now in order to realistically investigate the effect of an 
applied magnetic field on the two-dimensional system of 
antiferromagnetically correlated electrons we consider a 
case of half a flux quantum per plaquette, that is, ^ = 
2 = i. The energy dispersion relation is then obtained 
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Here the staggered magnetization (AF order) m is af- 
fected by the applied field. The self-consistent mean field 
equations for m and fi are given by 
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where the prime denotes integration over the reduced 
Brillouin zone of (k x ,k y ) with — \ < k x , k y < 5, as shown 
in Fig. |(a). 
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FIG. 2. Phase diagram of the staggered magnetization 
(AF order) in the (T, S) plane with (a) U = At (b) U = 8t. 
The dashed and solid lines indicate the boundary of antifer- 
romagnetic states for the magnetic flux quanta per plaquette, 
2=0 and | = i respectively. Here the elementary flux 
quantum c/>o is set to be a unity. 
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In order to examine the dependence of the external 
magnetic field (the solid lines in Fig. ||) on the staggered 
magnetization (AF order) m in the T-S plane we solved 
numerically Eqs. (10) and (0). We find that the reen- 
trant behavior of the staggered magnetization (AF order) 
persists at a certain region of doping rate 6 despite the 
application of the external magnetic field. The reentrant 
behavior appears at smaller doping rates for weaker cor- 
relations (see the case of U — 4t), as shown in Fig. ^. 
The domain of antiferromagnetic phase is smaller com- 
pared to the case of zero external field (the dotted lines 
in Fig. ^). The band gap mU owing to the antiferro- 
magnetic correlations gets smaller in the presence of the 
external magnetic field as a consequence of reduced mag- 
netization m. Such effect of reduction was predicted to 
be much larger with U — At than with U = 8t, as shown 
in the figure. 

We now investigate the reentrant behavior of staggered 
magnetization more in detail. In Fig. ^|(a) we display the 
staggered magnetization m as a function of both temper- 
ature T and doping rate S. At exactly half filling, i.e., 
5 = 0, m reaches the maximum at zero temperature. We 
find that the temperature at which it reaches a maxi- 
mum value quickly increases and converges to a finite 
temperature near half filling, as shown in Fig. [|(b). The 
reentrant behavior with U — 8t occurs more clearly at 
a region of higher doping rates (say 0.42 < 5 < 0.48, as 
shown in Fig. ||(b)). We see more readily from Fig. ||(c) 
the variation of reentrant behavior with the increasing 
doping rate. We find explicitly the rapid decrease of the 
AF order (staggered magnetization) with doping rate at 
a given temperature, e.g., kT = li, as shown in Fig. ||(c). 

Let us now consider the limiting case of vanishing ex- 
ternal magnetic field (B = 0). We obtain from |5|) the 
following expression, 
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as a result of g = 0. The dispersion relation is then, by 
diagonalizing the Hamiltonian above Q , 
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The self-consistent equations for m and /i are, in form, 
the same as (|To| ) and ([y]) above. However the integral 
is now over the reduced Brillouin zone of (k x ,k y ) with 
\k x \ + \k y \ < 7r as shown in Fig. 0(b). 
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FIG. 3. (a) Three dimensional plot of the staggered mag- 
netization as a function of temperature T and doping rate S 
for U — 8t in the absence of magnetic field, (b) The solid line 
represents a boundary of the antiferromagnetic state. The 
dashed line denotes the temperature at which the staggered 
magnetization m reaches the maximum value as a function of 
doping rate 8. (c) Temperature dependence of the staggered 
magnetization (AF order) for various values of doping rate 5. 
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ally lower. Now for the limiting case of the zero external preprint, pond- mat /9602027 . 

field (the dotted lines in Fig. ||), we find that there exists 
a reentrant behavior with a larger domain of staggered 
magnetization (AF order) in the T-S plane compared to 
the case of the applied external field (the solid lines in 
Fig. |^). In qualitative agreement with our present re- 
sults, Halvorsen et al. @ also showed the reentrant be- 
havior using a different approach. They used the Hub- 
bard model within the self-consistent second-order weak 
[/-perturbation treatment, thus allowing to study only a 
limited region of small U. Lately, Inaba et al. JtJ found 
qualitatively a similar reentrant behavior using the t-J 
Hamiltonian in the slave-boson representation, permit- 
ting to investigate only the large U limit case. On the 
other hand, our present approach allows to examine the 
full range of U, although we are not able to display all 
details due to a limited space. Further, unlike our present 
case, their studies were limited only to the case of zero 
external magnetic field. 

In the present study we derived the generalized 
Harper's equation which incorporates correlation effects 
between electrons. Unlike other studies [|]|7j we were able 
to examine as a function of both temperature and dop- 
ing rate the variation of the staggered magnetization m 
for the system of antiferromagnetically correlated elec- 
trons under the external magnetic field. By deducing the 
phase diagram in the T-S plane, we explored staggered 
magnetization m as a function of both temperature and 
doping rate and found the reentrant behavior even in the 
presence of the applied magnetic field. Finally, it is of 
note that the staggered magnetization or the AF order 
for other values of magnetic flux can be investigated 
in a similar manner. Further studies regarding this prob- 
lem will be resumed on. 
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